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We study theoretically a bistable systems with time-delayed feedback driven by weak periodic 
force. The effective potential function and the steady-state probability density are derived. The 
delay time and the strength of its feedback can change the shapes of the potential wells. In the 
0^ ' adiabatic approximation, the signal-to-noise ratio (SNR) of the system with a weak periodic force 

is obtained. The time-delayed feedback modulates the magnitude of SNR by changing the shape of 
the potential and the effective strength of signal. The maximum of SNR decreases with increasing 
the feedback intensity e. When e is negative (or positive), the time delay can suppress (or promote) 
<—j ' the stochastic resonance phenomenon. 
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In recent years, many complex systems with time-delay feedback have been investigated in optic laser systems, 
neural networks, coupled oscillators biological control, economic market, etc. In many cases, the delay reflects the 
transmission times related to the transport of matter, energy, and information through the system under investigation. 
^ , Therefore, time-delay systems can be regarded as simplified but very useful descriptions of systems involving a reaction 
chain or a transport process. Noise and delayed time are two important elements in many complex systems. Such 
systems can be viewed as a special case of stochastic systems with memory, which have recently studied numerically [l], 
0] and analytically 0, 0, [EH- 

Stochastic resonance (SR) behavior is one of the most studied and utilized fundamental physical phenomena and 
^ ■ has already observed through experimental [7[ ■ The original work on SR is mentioned by Benzi et al. [8| , for explaining 
the periodic recurrences of the earth's ice ages. The vast majority of studies on SR focus on the two-state model for 
symmetric systems 0, [1(1 HH [Hi m the linear response theory (LRT). For a normal bistable system, the dynamics 
can be divided into two different regimes, the linear response regime {A/ D <C 1) and the nonlinear regime {A/ ' D 3> 1) 
(where A and D are the amplitudes of the external periodic field and the noise strength, respectively) . This point 
was discussed in d, (TJ, [lj| • Shneidman et al. [lj| studied the weak noise limit of SR in a bistable system. The SR of 
asymmetry systems have been also investigated widely by the authors of res. [ill 53, EM Ell • Later, Nikitin et al. [T3 | 
I O I • expanded a comprehensive theory to study the effect of asymmetry on the switching dynamics in a bistable system 

with the limit of nonlinear response (weak noise limit). 
£N) • Recently, the time-delayed bistable system with small noise and magnitude of the feedback, has been widely studied 
^ ' in theory and experiment, especially the investigations on the residence time distribution and power spectra, having 
a great progress [2(1 [U [22j]. In the case of large delayed time (r rs 250), Tsimring and Pikovsky calculated the 
correlation function and the power spectrum to predict the coherence resonance (CR) and SR behaviors [23|. The 
phenomenon of CR was verified by experimentation (24[. In addition, the studies on delayed noise (noise recycling) 
have also attracted more and more attention and already had great results [25l . [26j . 

In this paper, we study the model with small delay and weak periodic signal. This stochastic system is a delay- 
induced non-Markovian process. So it is hard to obtain a appropriate analytical result. Using the small-delay 
approximation of the probability density 0, [H, [27|, we derive the time-delayed Fokker-Planck equation and the 
effective Langvein equation. It is found that the SR indeed could be effected by delays from the analytical results for 
the effective potential function, the steady-state probability density, and the signal-to-noise ratio. 
A typical time-delayed bistable systems is (23[ 



^± = x (t)-x 3 (t) + ex(t-T) + V2Dm- (1) 

Here r is the delay and e is the strength of the feedback, which also can be regarded as a modulation parameter of 
the system itself. is the Gaussian white noise with (£(£)) = and (£(t)£(i')) = S(t — t'). D is the strength of the 
noise. 

We consider a classic model driven by a weak periodic signal with the small delay 



dt 



= x(t)- x A (t) + ex(t - t) + A cos(Qt) + V2D£(t) , (2) 
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where A is the strength of the signal (in general, we adopt A = 0.02). For the low-frequency signal, the response of 
the systems to the weak periodic force is approximately linear, and the adiabatic limit is reasonable approach. The 
presence of small delay makes the bistable potential wells shack, but it is not enough to cause coherent resonance. 

The dynamics in Eq. ^ is a non-Mar kovian process. Using the probability density approach, the non-Markov 
process can be reduced to a Markov process and the approximate time-delayed Fokker-Planck equation is [H 

dP(x,t) _ d[h eff P(x,t)} ^ d 2 P(x,t) 

at ~ dx + d* x • [6) 

Here the conditional average drift h e f / reads 

f b 

h e ff = / dx T h(x,x T )P(x T ,t — t | x,t), (4) 

J a 

where x T — x(t — r), h{x, x T ) — x — x 3 + ex T + Acos fit, h{x) — x — x 3 + ex + A cos fit. The integral boundary a, b tend 
infinite (±oo). P(x r ,t — t | x,t) is the zeroth order approximate Markovian transition probability density [El [27|. 

„, 1 / (x-r — x — htxW) 2 \ , . 

P(x T ,t ~r\ X ,t) = 7 = exp (- iDr [)) J . (5) 

Substituting Eq. © into Eq. g]), we obtain 

h eff = (1 + er)(x - x 3 ) + e(l + er)x + (1 + er)A cos(Oi). (6) 
So, the effective Langevin equation for Eq. ([3]) becomes 

= (x - x 3 ) + ex + A cos(ffi) + ct[(x - x 3 ) 

+ex + Acos(Q,t)]+V2D£(t). (7) 

Note that a coupling term er [(x — x 3 ) + ex + Acos(ilt)\ is produced for the presence of time-delayed feedback. It 
indicates that the system is modulated by the delay and its feedback. From the effective Langevin equation obtained, 
one can investigate many dynamical properties of this system. 

The unstable point xo and two stable points x± can be calculated easily. In the absence of the external periodic 
force, x± = \/l + e. It means that the delay cannot affect the positions of the unstable points and stable points. 

The effective time-delayed potential function of Eq. ([7]) is 

U eff (x) = -(1 + er) Qx 2 - ^x 4 ^ ~ ^e(l + er)x 2 + (1 + e T )Acos{^t)x, (8) 
and the steady-state probability distribution function (PDF) is 

P st =Nexp(- l ^Y (9) 

where iV is the normalization constant, N — (J exp^^efce) . 

In this model, the modulation parameter e is expected to be larger than —1 (or should be larger) to ensure the 
bistability of the time-dependent potential U(x, t, r) possessing minima, respectively, at x±. Moreover, the depth and 
width of the potential wells increase with increasing of the strength of e, since the particles need larger noise D to 
cross the barrier. In order to limit the size of D, e is also not too large. In general, the positive e is not more than 1 
in this paper. All viewpoints are discussed on this basis. 

Without the periodic force, we depict the forms of potential for different time delays and their feedbacks in Fig.QJa, 
b). It is found that the depth of the potential wells increases with increasing of time delays for e > 0, but the width 
between two wells (two peaks) remains unchanged. That is to say the peak of the PDFs becomes sharper for deeper 
potential wells by increasing the ability of capturing particles. It means that the peak values of P s t increase with 
increasing of time delays. However, for e < 0, the results are opposite. In Fig.[2{a, b), when e = 0.6, r = 0.1 and 0.3, 
the peak value P max from the analysis is (0.508, 0.539), and from the simulation is (0.505, 0.523). Correspondingly, 
for e = —0.4, P max is: from the analysis (0.684, 0.663), and from the simulation (0.681, 0.657). The results are 
consistent with the above discussion. Fig. QJb) shows that the depth and width of the potential wells increase with 
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FIG. 1: (a) Potential U e ff(x) for different delayed times with e = 0.6. (b) Potential U e ff(x) for different strengths of time-delay 
feedback with r = 0. All the plots remove the effect of the external periodic forces. 



increasing of the strength of time-delayed feedback e, therefore the particles need larger noise to cross the barrier (see 
Fig. 02). 

Fig. [5] and Fig. [3] show the steady-state PDFs from the analysis and numerical simulations and their standard 
deviation (SD) a. It is found that the small-delay approximation of the probability density is a good approximation 
in the case of small delay (corresponding to a < 0.01). For the smaller noise, the availability of the approach is much 
better. In Fig. [2j the steady-state probability density of Eq. ([2]) (the circles) is consistent with Eq. (|7|) (the solid line) 
for r < 0.4, especially within r = 0.3. 

In general, the dynamics of a bistable system can be divided into two different regimes, the linear response regime 
and the the nonlinear regime. In the limit A <^ D, the response of the system to the periodic force is approximately 
linear (regarding A as perturbation term) and the linear response theory can be applied. In our systems, both the 
transition time of the probability from one well to another one and the signal period are much longer than the 
relaxation time of the system. As a result, the adiabatic limit works valid too. 

From Eq. ([7]), note that the unstable and stable points of the system are not influenced by the time delay r. It 
indicates that the particle spends most of the time near the minimal potential x — ±Vl + e, occasionally jumping 
from one to another because of the noise. The transition rates between the wells are Kramers' escape rates |9|, [28| 



(l + £ )(l + er) 
7± = TTn ^ 



4D 



x exp 



±(1 + er)y/l + eAcosttt 



D 



(10) 



Here, 7+ is the transition rate of the probability from left-hand well to the right-hand, and vice versa. 

Normally, the periodic force is small (Ax <C AC/) and slow (CI <C loq), where is the frequency of system 
vibrations inside the well. In this paper, the frequency of signal Q, = 0.003. Without the periodic force, we have the 
approximation of Eq. HI 



(l + e)(l + er) 

7±0 = — 7T, — exp 



(l + e)\l + er) 



AD 



(11) 



Consequently, the Kramers' escape time is 7 1 . In general, the time delay r <C 7 1 - 

The power spectrum density of the output variables, a Fourier transform of the autocorrelation function, in Eq. ([7]) 
is given by [lllS EH 



where 



AT _l V2(l + e)(l + er) 
N = 7+0 + 7-0 = exp 



M' 2 



2(l + e)N 



2(N 2 + n 2 ) 

(l + ef(l + er) 
AD 



(12) 



M 



A 
D 



(1 + er)VT + eiV. 



(13) 
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FIG. 2: The steady-state probability density P st for different delays r with parameters D and e: (a), (b), (c) e = 0.6, D = 0.4; 
(d), (e), (f) e = —0.4, D = 0.08. The cycles represent the results obtained from simulations Eq. and the solid lines results 
form Eq. ©, A=0.02. 



In Eq. (HU), < 1 allow us to omit the last term -Af 2 / [2(N 2 + fi 2 )] . So S(w) is denned only for positive fi. Then, 
Eq. IT2|) becomes 

where S±(u>) is the output power spectrum of the signal, and ^(w) is of noise. From the definition of SNR, i? = 
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FIG. 3: The standard deviation (SD) cr = ^/TJJn — 1) X)r=i(P ia ~~ Pis) 2 - P»a i s the analytical value of P s t at the ith point, and 
p i3 is the average value of P 3t from simulation at ith point. 



Si (ui) / S2 (u> ) , we have 



^= + 6)2(1 + £T)3eXp 



(l + £ ) 2 (l + e r) 



4L> 



(15) 



From dR/dD = 0, we get 



D 



;(1 



(1 + er) 



-AC/. 



(16) 



Here, AU is the height of potential barrier. When D/AU « 0.5, the SNR reaches a maximum. 

In Fig. [H we show the curves of SNR as a function of the noise intensity D for different e and r. Note that each 
curve exhibits an optimum noise intensity where the SNR has a maximum, which is the characteristic signature of SR 
phenomenon. Furthermore, it is obvious that the peak value of SNR decreases with increasing of feedback e when r 
is a constant [see Fig. Ufa)]. In the case of negative feedback, the time delay suppresses the SR phenomenon, on the 
contrary, the larger delay will help to enhance the peak of SNR for the positive feedback [as shown in Fig. (Hb)-(c)]. 

It is interesting to compare the effects of delay and its feedback on R max (the maximum of SNR ). When D/AU = 
0.5, the maximum of SNR is 



Rn 



16V2e~ 2 A 2 



(1 + e) 2 



(17) 



The model [Eq. @] is controlled within a small delay. Note that e = is the critical point. The negative e is more 
effective to promote the SR than the positive [see Fig. [UJa)]. In Eq. (fT7|) and Fig. [5{b), the linear relationship of 
Rmax and r is displayed when e is fixed. Obviously, within a small delay, it is appropriate to observe the clear SR 
phenomenon when the feedback is negative. 

Although the delay r <C 7" 1 and r <C T (where 7 _1 and T are Kramers' escape time and the period of the signal), 
the effect of time delay on SR cannot be neglected. Eq. ([2]) can be rewritten by the following form 



dx(t) 
dt 



(1 + e)x(t) - x 3 {t) + e[x(t - r) - x(t)] + Acos(Slt) + V2D£(t), 



where e[x(t — t) — x(t)] is a small quantity. We define 

e[x(t - t) - x{t)\ = e X (r)x(t), 



(18) 



(19) 



as the function of r. x( T ) 1S a small quantity too, and is a linear relationship with the delay. Then the another form 
of Eq. d2) is presented, 



dx(t) 
dt 



[l + e + ex(r)]x(t) - x 3 (t) + Acos(ftt) + V2D£(t). 



(20) 
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FIG. 4: SNR as a function of the noise intensity D/AU with the magnitude of the signal A — 0.02, where AU is the height of 
the potential barrier without external periodic force. Parameters are: (a) r — 0, (b) e = —0.4, (c) e = 0.6. 
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After the moment of a particle crossing the barrier, from the left well to the right, x(t — r) < x(t), x(t) > 0, then 
x(t) < 0; contrarily, when the particle just jump into the left well, x(t) < 0, x(t — r) — x(t) > 0, and %(r) < 0. 

It is well known that the depth of the potential wells AU (or the potential barrier) is one of the important parts 
to contribute to SNR. In our model, the amplitude of the signal A is a constant, thus its ability to modulate wells 
will keep unchanged. When AU increases with e, the particles are difficult to cross the barrier, which prevents the 
noise and signal from achieving resonance. Besides, with the enhancement of the noise background S*2(w) (for keeping 
D I AU 0.5), SNR will weaken. It means that the peak of SNR will decrease with the raising of e, which is shown 
in Fig. gK a). 

When the delay r exists in this system, at the moment of a particle just crossing the barrier, the function x( T ) < 0. 
For the negative feedback e, the term ex{r) > 0, 1 + e + e%(r) > 1 + e, which causes the potential barrier AU larger 
in Eq. (p?0"|) , thus the peak value of SNR will smaller [see Fig. 0|b)] . That is to say that the delay r prevents the 
particle returning to the former well after it just reach one well, and this behaviors lead to SR phenomenon weakening. 
Whereas the term ex(r) < with the positive e, 1 + e + ex( r ) < 1 + e, the potential barrier AU is depressed, and the 
particle can easily return to the original potential well. So r is helpful to raise the peak of SNR [shown in Fig. S(c)]. 

From the part induced by the time delay r in Eq. ([7]), er [(x — x 3 ) + ex + Acos(Qt)] , the external periodic signal 
and the original system are modulated together by the time delay r and its feedback intensity e. Although r can tune 
the potential barrier AU from the negative (or positive) feedback e, the amplitude of the signal A* is also modulated 
by e and t, A* = (1 + er)A. So the variations of SNR displayed by Fig. [1] are proper. 

Generally speaking, in the stochastic system driven by weak periodic force, because of the presence of the delay t 
and the strength of its feedback e, the system produces a coupling term, such as er [(x — x 3 ) + ex + Acos(fii)] in this 
work. By changing r and e, the signal amplitude A is modulated, and the peak value R m ax of SNR exhibits different 
magnitudes under the same intensity of noise (D/ AU « 0.5) and shows different tendency. When r keeps unchanged, 
the SR phenomenon of the negative e (—1 < e < 0) is more easily observed than it for the positive e. For the positive 
e, the SNR of the system is promoted with increasing of the time delay r. On the contrary, the SR phenomenon with 
the negative e is suppressed by r. 

In this paper we study a simple time-delayed bistable model with Gaussian white noise driven by weak periodic 
force. Using the small delay approximation of the probability density function, the time-delayed Langevin equation 
is extended to a effective Langevin equation. Correspondingly, we obtain the effective potential function U e ff and 
the steady-state probability density p s t . For a weak periodic force, the analytical expressions of SNR and R m ax are 
derived in the adiabatic limit, and the influences of delays on them are also discussed. It shows that the time delay 
is one of the most important elements to affect the dynamics of complex systems. 
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